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Detection of Two Massive CO Systems in 4C 41.17 at z = 3.8 
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Abstract. We have detected CO(4-3) in the z=3.8 radio galaxy 4C 41.17 with the IRAM Interferometer. The CO is in two 
massive (M(jy„ ~ 6xl0"'Mo) systems separated by I'.'S (13kpc), andby 400kms"' in velocity, which coincide with two different 
dark lanes in a deep Lya image. One CO component coincides with the cm-radio core of the radio galaxy, and its redshift is 
close to that of the Hell A 1640 AGN line. The second CO component is near the base of a cone-shaped region southwest of 
the nucleus, which resembles the emission-line cones seen in nearby AGN and starburst galaxies. The characteristics of the 
CO sources and their mm/submm dust continuum are similar to those found in ultraluminous IR galaxies and in some high-z 
radio galaxies and quasars. The fact that 4C 41.17 contains two CO systems is further evidence for the role of mergers in the 
evolution of galaxies at high redshift. 
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1. Introduction 

_ Because of its luminosity and large angular extent, 4C 41.17 
. has become the most studied high redshift radio galaxy 
■ (HzRG). //primages show that its host galaxy contains several 
star-forming components, including (i) a linear radio-aligned 
feature wit h spectroscopic c haracteristics of a young stellar 
population jPev et all Il997 >). (ii) a clumpy system, separate 
from the radio source, and (iii) low surface brightnes s UV 
emission extending over 70 kpc 6van Breugel et allll999l). The 
entire system is embedded in a giant Lya halo /Reuland et all 
I2OO3I) and diffuse soft X-ray emission ( Sch arf et al., 2003). In 
spite of the detection o f its dust emission (IDunlopet all 1 19911 
IChin i&Kriige?.'l994^ searches for molecular gas have been 
unsuccessfu l up to now (Ivison et al., 1996; Evans et al., 1996; 
iBarvainis & Antonuccilll99dl:IScoville. Yun. & Brvantt ll996). 
In this letter, we report the first detection of two massive CO 
components near the centre of this large forming galaxy'. 



2. Observations and data reduction 

We observed 4C 41.17 at 3.1 and 1.2 mm with the IRAM 
Plateau de Bure Interferometer in its compact D configuration 
between 1999 August and 2003 August. The spectral correlator 



' We adopt Ho=7 1 km s"' Mpc"' , nM=0.27 and fiA=0.73. At z=3.8, 
the luminosity distance £'/.=34.4Gpc, and 1" corresponds to 7.2 kpc. 



covered a bandwidth of 580 MHz at each frequency. At 3. 1 mm, 
the 6-antenna equivalent on-source observing time was 29 h, 
and the beam was 7'.'0 x 4'.'7 at PA 107°. The central fre- 
quency was initially tuned to 96.090 GHz to observe CO(4-3) 
(vrest = 461.040 GHz, but later re-centred to 96.250 GHz to 
better cover the detected line. After phase and amplitude cal- 
ibration, we merged the data at different centre frequencies 
into a data-set of 36 channels of 20 MHz each. As a result, 
the outer 7 channels on each end have cra;0.9mJy compared 
to cr»;0.6mJy in the central ones. In the final data cube, the 
zero point of the velocity scale is 96.093 GHz, cor responding 
toz^ 3.79786, the redshift of the He II A 1640 Une jPev et all 
I1997I) . We made naturally weighted maps using the AIT'S task 
IMAGR, with CLEANing appUed only to the velocity-integrated 
maps. 

We observed simultaneously at 1.2 mm (241.453 GHz 
DSB) to study the dust continuum. We used only 1 .2 mm data 
taken with a precipitable water vapour content < 3 mm, which 
gave a usable on-source observing time of 7.7 h, and an rms 
noise of 0.8 mJy/beam. The 1.2 mm beam was 2'.'9 x r.'5 at PA 
100°, but we convolved the image with a 3'.'0x3'.'0 Gaussian. 

In 2004 March, we also observed the 4C 41.17 field at 
1.2 mm (-250 GHz) with the 1 1 7-element MPIfR M illimeter 
Bolometer array (MAMBO-2; Krevs a et all Il998l) at the 
IRAM 30m telescope. The beam FWHM is W7 with an array 
size of 4'. We made eight on-the-fly maps, with 41 subscans of 
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Fig. 1. CO (4-3) spectra at the positions of the red (top panel) 
and blue components (bottom). The thin green line shows 
the data in 62 kms ' channels, and the thick black line in 
184 kms"' channels. Error bars are Icr, and velocities are rel- 
ative to 96.093 GHz (zlsr = 3.79786). The red and blue bars 
indicate the linewidths; the thick part of the red bar marks chan- 
nels used in Fig. 3, left. Note that these two positions are less 
than a beam size away, and are therefore not fully independent. 



Table 1. Observed and derived parameters for CO(4-3) in 
4C 41.17. 



Parameter 


Red component 


Blue component 


Total 


RA(J2000)° (peak) 


6''50"'52'.24 


6''50"'52n2 


6''50'"52^17 


DEC(J2000)° (peak) 


41°30'3I'.'6 


4r30'30'.'3 


41°30'30'.'9 


/ Sco<^V[Jy kms-'] 


1.2±0.15 


0.6 ±0.15 


1.8 ±0.2 


Central velocity [ km s"' ] 


-130 ±50 


-550 ± 100 


-285±100 


Velocity width [ km s"' ] 


500 ± 100 


500 ± 150 


1000±150 


L^o(4-3) [Kkms^'pc^] 


44x10'° 


2.2x10'° 


6.7x10'° 


M(H2) [Mo] 


3.6x10'° 


1.8x10'° 


5.4x10'° 



" Positional uncertainty 0:03 in RA and ^.'3 in DEC 
>> Relative to 96.093 GHz (z =3.79786). 



40 s each, while chopping the secondary mirror in azimuth at 
2 Hz by 39, 4 2, or 45". We reduced the data using MOPSIC 
llZvlkalll998h . The map covers 3'x3' with an rms noise of 0.8 
mJy. 

3. Results 

Figure [2 shows CO(4-3) spectra at two positions. No line 
emission is detected in the range -1550 to -1000 kms"' and 
-1-250 to -1-800 kms"' to 0.75 mJy (3o-). In these outer chan- 
nels, we find a marginal continuum emission at ~0.3mJy, 
which is consistent with the llmJy at 850jum jlvison et all 
I2OOOI) and 3.8 mJy at 1.2 mm (our MAMBO map) extrapo- 
lated to 3.1mm (~0.2mJy), and the non-thermal contribu- 
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Fig. 2. CO (4-3) position-velocity slice at PA=51° (i.e. along 
the radio axis). Horizontal Unes mark the radio core and outer 
lobes. The position zero is 6''50"52n7 41°30'30'.'9, the beam 
along this cut is 6", and the contour steps are 0.4 mJy, with the 
first contour at 0.4 mJy. 



tion of <0.2mJy (over the 13" source), extrapolated from 
6.4 mJy at 2 cm (Ch ambers et al., 1990) and 2.7 mJy at 1.2 cm 
jlvison et all ll996'). Because of these low values (below the 
first contour in Fig.|5J, we do not correct the CO fluxes for the 
continuum. 

The position- velocity slice (Fig. O shows that the CO 
emission has two components: (i) a 'red' component at ~ 
-130km s"' relative to z=3. 79786, and (ii) a 'blue' compo- 
nent at -550 km s"' . Figure|3lshows the integrated red and blue 
components^, separated by 1"8 (13 kpc projected). Table^lists 
the observed parameters, the line luminosity L^q(4-3), and 
the molecular gas mass M(H2), calculated assuming a constant 
brightness temperature from CO(4-3) to CO(l-O) and a con- 
version factor Xco=0.8 Mo(K kms~'pc^)"' derived for lo cal 
ultraluminous infrared galaxies jPownes & Solomon. 1998|). 

Figure ^ shows three maps of the dust continuum in 
4C 41 . 17. Unlike the solid lOcr and 5cr detections in the 850 jum 
SCUBA and 1 .2 mm MAMBO maps, our 1 .2 mm interferom- 
eter map shows a 4.3cr peak at 06''50"5Z24, -h41°30'31'.'9 
(J2000; i.e. at the red CO component) only after convolution 
with a 3'.'0 Gaussian. This PdBI map yields a flux 5 1.2mm - 
3.0 + 0.7 mJy, consistent with the S 1.2mm = 3.8 + 0.6 mJy from 
the MAMBO map. The S/N in our 1 .2 mm maps is insufficient 
to constrain the spatial extent of the thermal dust emission re- 
ported by Ivison et al. (2000) and Steven s et al. (2P_03). 

Both CO components are gas-rich systems with Mhj ~ 
3 X 10'" Mq (TableQJ. Their projected separation of 13 kpc and 
apparent relative velocity of ~400 km s"' imply a combined dy- 

- To increase S/N, the red component in Fig. |3| covers only the 
brightest 250 km s"' , indicated by the thick bar in Fig.Q(top). 
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Fig. 3. Velocity-averaged CO maps (thick/white+black contours) of the red (left) and blue (right) components. CO contours are 
-3, -2, 2, 3, lOo", with cr=0.2 mJy beam The lower left inset shows the CO beam. The two CO maps are superposed on the 
Lya image, each with its own colour scheme to show the two dark lanes in the Lya. The thin green contours (left panel) show 
the 1.4 GHz radio map (Carilli et al. 1994), and the yellow cross indicates the radio core. The thin straight lines (right panel) 
indicate a possible AGN or starburst emission-line cone. 



namical mass RAV^/G of > 6 x IO'^Mq, an order of magnitude 
smaller tha n the typical baryonic mass of HzRGs integrated 
over 64 kpc jRocca-Volmerange et allEo04 l). 

We associate the red CO component with the AGN because: 
(i) it agr ees within the uncertainties, w ith the cm-radio core po- 
sition tCarilh. Owen. & Harris! fl994h and the hard X-ray point 
source ( lsclwFetaLr]2003l) . (ii) it coincides with the central 
dark lane seen in the Lya image (Fig. 01 left, near the cross) and 
the HST /?-band image (.Bicknell et aL .2000,) . and (iii) it has 
the same velocity offset as the low-ionizati on interstellar ab - 
sorption lines in the deep Keck spectrum o f bev et alJlll997[) . 
The CO redshift of 3.7958±0.0008 is thus likely to be the real 
systemic redshift of the host galaxy of 4C 41 . 17. 

The blue CO component peaks at the position of a second 
apparent dark lane in the Lya image (Fig.|3j right), which sug- 
gests that the Lya may be absorbed by dust of another kine- 
matic system, associated with this CO. This apparent Lyes' ab- 
sorption lane is in a region which resembles the emission-line 
cones seen, although on smaller scales (a few kpc), in some 
nearby AGN (e.g. Cygnus A; Canalizo et al., 2003) as well as 
in starburst galaxies (e.g. iVeilleux & Rupke...200Z) . 

4. Discussion and Conclusions 

What is the relation between the molecular CO-gas, the dust, 
the Lya halo, and the massive forming galaxy with its radio- 
loud AGN? Our CO velocity profile (Fig.EJl shows a remark- 
able similarity to the Lya velocity profile of .Dev.Cl999.) . Both 



the CO and Lyes' are split into two components, separated by a 
projected distance of ~13 kpc, and by ~400 kms"'. This sug- 
gests that the Lya emission may also come from two separate 
components (rather than a single component split by an asso- 
ciated HI absorber). Ho wever, Lya traces mu ch less dense gas 
(«H ~ 17-150 cm'^: .Villar-Martm et allE003 ) than CO. which 
must have a density of 10^ to lO'^cm"^ to have enough CO-line 
opacity to give a typical brightness temperature of order 30 K. 
Hence, the CO and Lya may originate from the same gas-rich 
systems, but they do not necessarily trace the total extent of 
these regions. This is obvious from the 200 kpc spatial extent 
of the Lya, while the CO emission is unresolved with our 6" 
(43 kpc) beam. If fact, we can put even stronger constraints on 
the size of the CO-line sources, using some basic assumptions. 

The FIR dust luminosity of 4C41.17 is very high, Lfir ~ 
1.5 X lO'-' L0, a value typical of an ultraluminous starburst. The 
FIR fluxes, including our new continuum data points at 3. 1 and 
1.2 mm, imply a dust temper ature of 54 ± lOK, in agr eemen t 
with eariier values (Benford eTail Il999t Hcharf et all l2003h . 
Although the mm/sub-mm continuum is optically thin, the CO 
lines are not, so the observed brightness temperature of the 
CO will be about the same as the gas temperature. The ex- 
isting interferometer CO maps of ULIRGs, the only nearby 
objects with comparable FIR luminosity, show that the bright- 
ness temperatures of the low-J CO lines are compara ble with 
the FIR dust temperature jPownes & Solom'onl ll998V Indeed, 
for high-z CO detections even to be possible, the gas must 
have a significant brightness temperature, typically 30 to 50 K, 
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Fig. 4. Thermal dust emission in 4C 41.17. The colourscale 
shows the 1.2 mm PdBI map convolved with a 3'.'0 Gaussian 
(not corrected for the 20" primary beam). Superposed are the 
850//m SCUBA map (thin black contours; beam 14") and 
our 1.2 mm MAMBO map (thick white contours; beam 11"). 
Contours are -2, 2, 3, 4, 5, 6, 7, 8, 9 and lOcr, with cr=1.2 and 
0.8 mJy beam"^ for the SCUBA and MAMBO maps, respec- 
tively. The white cross indicates the cm-radio core. 



over several hundred pc. For 4C41.17, this means that if the 
CO(4-3) brightness temperature is ~30K, then the observed 
CO(4-3) luminosity (Table 1) implies a CO source diameter 
d=(4 L'^Jin Tb AV)f-^ of 1 .4 to 1 .8 kpc, or ff.'l to 0'.'25. 

This size and H2 mass imply a hydrogen column density 
of order lO^'* cm"^, which is consistent with the observed mm- 
FIR dust spectrum becoming opaque near restframe 100//m. 
Applying the Stefan-Boltzmann formula with the derived CO 
source diameter ii~1.6kpc and dust temperature then yields 
■^'FiR - n<T d^T'^^lO^^ Lq, similar to the observed value. The 
CO source size is also roughly consistent with the total mass 
and typical density of the CO gas (10^ to 10'*cm"^), which im- 
ply source sizes of order ~ 1 kpc, depending on geometry. Note 
that these are maximum sizes of the CO. If one assumes the CO 
is in more than two sources, then each component will have a 
smaller size. We also note that our observations are insensitive 
to more widely distributed, cooler CO. 

The only known place where such high gas densities over 
such dimensions are found are the circumnuclear disks ob- 
served in ULIRGs and some quasars. The 'red' CO source 
in 4C41.17 thus finds a natural interpretation as a circumnu- 
clear starforming disk around the radio-loud AGN. The fact 
that there appears to be a second such CO source, 13 kpc away, 
separated by ~400 km s"' in velocity, that also contains about 
the same mass of molecular gas, suggests the AGN and possi- 
bly starburst activity may have been triggered by the interac- 
tion of the two objects. The absence of UV/optical continuum 
signatures of starburst activity at the position of the blue CO 



component remains surprising. A possible explanation could 
be that th is starburst has n ot yet reached its peak UV/optical 
emission faaas et alll2003l) . 

To summarize: our observations indicate that the 4C 41.17 
system contains two massive CO components, each of which 
may be associated with an obscured black hole. This is 
remar kably similar t o the t wo CO systems in 4C 60.07 
( Papa dopoulos et all I2OO0I: iGreve. Ivison. & Papadopoulol 
i2004l) . There are three other HzRGs, 6C 1909-1-72, 
B3 J2330-H3927 and TN J0120-H1320 w i th detec ted CO 
jPapadopoulos et all l200dl lOe Breuck et all l2003albl) . It will 
be of interest to determine if they are also double sources 
in CO, which would further indicate the role of mergers in 
triggering AGN activity in the most massive galaxies at high 
redshift. 
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